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1.0  INTRODUCTION 


This  report  presents  an  improved  theory  for  the  analysis  of  input  impedance  of 
a  probe-fed  microstrip  patch  antenna.  The  motivation  of  this  work  arises  from 
the  need  for  a  more  accurate  determination  of  the  characteristics  of  millimeter 
microstrip  antennas  than  obtainable  from  currently  available  cavity  and  trans¬ 
mission  line  models.  A  complete  spectrum  representation  of  a  probe  excited  micro¬ 
strip  antenna  has  been  developed  based  on  the  work  of  Chang  and  his  colleagues 
[1,3],  where  the  cavity  wall  reflection  coefficient  has  been  obtained  as  a  func¬ 
tion  of  the  angle  of  incidence  of  interior  waves  bouncing  between  the  walls. 

The  first  mathematical  model  of  the  basic  microstrip  radiator  used  a  one-dimen¬ 
sional  transmission  line  analogy  [4,5]  and  were  later  extended  to  two  dimensions 
by  the  cavity  model  [6-8].  However,  neither  of  these  theories  account  for  the 
dispersion  effects  caused  by  the  grounded  dielectric  slab.  Of  particular  interest 
here  is  the  background  work  of  the  modal  field  expansion  theory  [6],  This  theory 
is  based  on  the  conventional  Fourier  series  expansion  of  the  fields  inside  the 
cavity  formed  by  the  microstrip  patch  and  its  associated  ground  plane.  In  its 
present  form  it  assumes  a  thin  cavity  so  that  the  interior  fields  are  composed 
of  discrete  modes  which  may  be  calculated  assuming  a  fixed  angularly  independent 
wall  reflection  coefficient.  The  effect  of  radiation  is  represented  in  terms  of 
an  increased  substrate  loss  tangent  as  done  by  Lo  et. al.  [6]  or  by  the  method  of 
impedance  boundary  condition  at  the  wall  by  Carver  et.  al.  [8].  When  the  electrical 
thickness  of  the  cavity  ranges  from  0.02  to  0.2  as  in  the  case  of  practical  micro¬ 
strip  antennas  in  the  L,  S  or  C  band  regions,  the  surface  waves  supported  by  the 
dielectric  substrate  carry  relatively  little  power  and  the  assumption  of  a  per¬ 
fect  magnetic  wall  and  an  angularly  independent  reflection  coefficient  for  any 
given  patch  geometry  and  frequency  of  operation  is  fairly  reasonable. 

By  contrast,  microstrip  patch  antennas  constructed  with  commercially  available 
substrate  thicknesses  and  operated  at  the  X-band  or  higher  frequencies  may  have 
an  electrical  thickness  equal  to  and  exceeding  0.5.  Experimental  results 
obtained  with  the  X-band  patch  antennas  [9]  show  that  the  cavity  theory  breaks 
down  as  the  substrate  becomes  electrically  thicker  since  an  appreciable  amount 
of  input  power  is  transferred  to  surface  waves  on  the  open  grounded  substrate. 


This  inadequacy  of  the  cavity  theory  underscores  the  need  for  a  better 
theoretical  description  of  millimeter  wave  microstrip  antennas. 

The  analysis  of  the  microstrip  patch  antenna  by  Change,  Chang  and  Kuester 
[1,2]  is  based  on  the  recognition  that  the  natural  modes  of  the  antenna  are 
established  on  transverse  resonance  conditions  which  incorporate  the  angularly 
dependent  reflection  coefficient  associated  with  the  cavity  walls.  This 
angular  dependence  is  a  direct  result  of  the  surface  waves  and  radiation  sup¬ 
ported  by  the  grounded  dielectric  substrate.  The  Wiener-Hoff  technique  has 
been  used  to  find  the  reflection  coefficient  and  the  analyses  has  been  used 
to  compute  the  optimum  range  of  patch  side  length  aspect  ratios  in  the  sense 
of  low  Q  or  greatest  bandwidth. 

The  present  theory  considers  the  complete  spectrum  of  plane  waves  propagating 
from  the  excitation  probe  at  all  angles.  The  natural  modes  are  established  on 
the  same  transverse  resonance  condition  obtained  by  Chang  which  incorporates 
the  angularly  dependent  reflection  co-efficient  associated  with  the  patch 
boundaries.  The  reactive  part  of  the  probe  impedance  for  a  truncated  rectan¬ 
gular  patch  is  due  to  the  evanescent  waves  which  are  usually  confined  to  the 
vicinity  of  the  probe.  Those  directed  at  specific  angles  interfere  construc¬ 
tively  and  propagate  down  the  strip  without  decay.  These  are  the  guiding 
modes  and  the  corresponding  edges  of  the  patch  act  as  the  non-radiating  walls. 
Some  of  the  spectrum  propagation  will  leak  out  at  the  edges  of  the  patch  in 
the  form  of  radiation. 

This  report  presents  the  theoretical  analysis  of  input  impedance  to  probe  fed 
microstrip  patch  antennas.  Expressions  have  been  obtained  as  a  function  of 
patch  size,  probe  location  and  size  and  frequency  of  operation. 


2.0  INFINITELY  EXTENDING  MICROSTRIP  PATCH  ANTENNA 


2. 1  Geometry  of  the  Problem 

The  analysis  has  been  broken  into  two  phases.  In  the  first  phase,  the  micro¬ 
strip  patch  antenna  has  been  considered  as  infinitely  extending  in  the  x 
direction  as  shown  in  Figure  1.  The  width  of  the  patch  is  i.  in  the  y  direction. 
It  is  excited  by  a  z-directed  probe  located  at  (x0,  yQ) .  The  patch  is  located 
on  the  surface  of  a  grounded  dielectric  slab  of  thickness  d  and  relative 
permittivity  e^.  This  then  allows  us  to  view  the  structure  as  an  infinitely 
extending  waveguide  and  helps  us  to  establish  all  waves  both  propagating  and 
evanescent  which  emanate  from  the  exciting  probe  at  different  angles. 

The  evanescent  waves  directed  at  different  angles  would  decay  very  quickly 
within  a  small  distance  from  the  probe  and  would  account  for  the  reactive  part 
of  the  input  impedance  (Figure  2a) .  Some  of  these  waves  may  reach  out  to  the 
walls  in  the  x  direction  where  they  could  partially  leak  out  as  spurious  radia¬ 
tion  and  also  be  partially  reflected  in  the  y  direction. 

The  propagating  modes  bounce  back  and  forth  in  the  y  direction  (Figure  2b) , 
either  interfering  constructively  or  destructively.  Only  those  directed  at 
specific  angles  interfere  constructively  and  propagate  down  the  patch  as 
guiding  modes.  They  ultimately  characterize  the  resonant  properties  of  the 
patch.  The  others  which  interfere  destructively  ljak  out  as  spurious  radia¬ 
tion  in  the  y  direction  as  the  waves  progress  along  the  x  direction. 

Hence  an  analysis  of  an  infinite  patch  would  help  establish  the  transverse 
resonance  condition  for  the  guided  modes  and  the  entire  spectrum  of  plane 
waves  emanating  from  the  source.  In  the  second  phase  of  the  analysis,  we 
consider  a  truncated  patch  which  results  in  waves  bouncing  in  tne  x  direction. 
The  input  impedance  of  such  a  finite  patch  can  then  be  calculated  from  the 
electric  field. 

2.2  Field  Under  the  Infinite  Microstrip  Patch 

We  first  assume  the  electric  current  density  due  to  the  filamentary  exciting 
probe  source  is  of  the  form 
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Fig.  1.  Infinitely  long  probe-fed  microstrip  patch  antenna. 
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The  primary  z-directed  electric  field  emanating  from  the  probe  then  satisfies 
the  wave  equation 
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where  n2  =  and  k2  =  u2p0e0  and  where  the  assumed  but  suppressed 
time  variation  is  exp  (-io-t).  It  is  well-known  that  the  solution  of  eq.  (2) 
can  be  represented  by  a  spectrum  of  plane  waves  of  the  general  form: 
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A  specific  solution  of  eq.  (2)  is 
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Eq.  (4)  represents  an  angular  spectrum  of  plane  waves,  or  primary  waves  from 
the  probe.  The  exponential  term  in  the  integral  (4)  can  be  written  for  prop¬ 
agating  (non-evanescent)  waves  as 
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where  r  is  the  radial  distance  from  the  probe  and  a  =  n  cos  $  with  the  angle 
<)>  defined  as  shown  in  Figure  26.  However,  a  is  in  general  complex  since 
evanescent  waves  are  also  possible. 

Reflection  occurs  as  these  primary  waves  are  incident  on  the  two  edges  of  a 
microstrip  of  finite  width.  For  an  infinitely  long  patch  the  waves  travelling 
in  the  x  direction  while  bouncing  in  the  y  direction  would  be  the  sum  of  inci¬ 
dent  waves,  those  which  travel  in  the  positive  y  direction  and  those  which 
travel  in  the  negative  y  direction.  Hence  we  can  write 
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The  amplitude  A  of  waves  reflected  from  the  edge  at  y  =  £/ 2  can  be  determined 
by  the  boundary  condition  imposed  at  the  edge: 

A  =  y0  +  B  £/zj  r(a)  eikjn*-a*  £/ 2 

where  r(a)  is  the  complex  reflection  coefficient  obtained  in  [2]. 
the  amplitude  B  is  given  by  the  edge  condition  at  y  =  -£/2  as 

B  =  [>»*  _  y0  +  a  eik^  -  a*  ifl\  r(a)  giWn2  -  a*  £/ 2  (8) 
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Similarly, 
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Solving  for  A  and  B,  the  sum  of  the  second  and  third  terms  in  eq.  (6)  would 
be  given  as  follows 
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Substituting  in  eq.  (6)  we  have 
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1 

l 

I 
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2.3  Angularly  Dependent  Reflection  Coefficient 


As  mentioned  earlier,  T(a)  is  the  complex  reflection  coefficient  of  the 
plane  waves  bouncing  in  the  y  direction.  When  the  plane  waves  from  the 
exciting  probe  proceed  in  the  y  direction  at  different  angles,  whether 
these  waves  radiate  into  open  space  or  not  depends  on  the  angle  of  incidence. 

A  complete  reflection  is  possible  if  the  angle  of  incidence  <f>  =  sin  ^  (a/n) 
where  n  is  the  refractive  index,  is  greater  than  some  critical  angle.  Beyond 
this  critical  angle  the  reflection  coefficient  has  a  magnitude  of  unity. 

The  reflection  coefficient  was  obtained  by  Chang  and  Kuester  [2]  using  the 
Weiner-Hopf  technique  as  applied  to  coupled  integral  equations  for  charge 
and  longitudinal  distributions  on  the  patch.  An  analytical  expression  obtained 
by  them  [3]  for  thin  slabs  is  given  by 

r(ct)  =  e“ix(ot)  (13) 


where 


X(a)  *  2  tan  1  (  -  -  tanh(A)J  -  f  (-/i2-a2) 

V  ^2^?  )  e 

A  S  ["teCv'SM'  kd)  +  y  -  lj 


+  2  v-v  ~  2  w 


fe(-^n2-a2)  5  ~2— n-2Jp2  -M  Jj_  jin(»^2-l'  kd)  +  y  -  lj 


+  2  Qo(-5e)  -  ^n(27T) 


(14a) 

(14b) 


(14c) 
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When  a~n,  A  is  very  small  and 


(l-a2)pr  |  £n(/a2-i  kd  )  +  ylj 


+  n2  |2Qo(-6e)  -  S,n(2ir)J 
-  a2  |2Qq(  <5p  )  -  £n(2TT)j 


e  -  1  Vi  -  1 

6e  =  — -  ;  6p  = — - -  ;  Q  (z)  =  £  zm£n(m) 

^  i-i  ..11  o  m—  X 


e  +  1 
r 


\  +  1 


(14d) 

(15) 


and  y  =  Euler's  constant  0.5772.  Equation  (14)  is  valid  for  kd  <  1  as  is 
found  in  many  practical  cases.  For  thicker  slabs,  one  can  employ  the  integral 
representation  given  in  [2]  directly. 


Physical  Interpretation  and  Analysis  of  The  Field  Expression 


The  transverse  resonance  condition  for  a  conducting  strip  of  finite  width  l 
occurs  when  the  total  phase  change  for  the  TEM  wave  bouncing  back  and  forth 
between  the  two  ends  of  the  strip  -2x(<*)  +  2k/n2-oi2  Z  is  equal  to  an  integer 
multiple  of  2tt ;  that  is. 


-2x(or)  +  2  kSjn^-a2  =  -2p n 


(16a) 


The  first  term  is  the  phase  accumulation  resulting  from  the  complete  reflection 
of  this  wave  from  the  two  edges  and  the  second  term  is  the  phase  change  across 
the  transverse  y  dimension  for  the  zig-zagging  wave. 


The  particular  values  of  a  =  a^,  p  =  1,  1,  2,...  that  produce  this  resonance 
correspond  to  the  modes  on  a  microstrip  structure.  These  can  be  determined  by 


a  numerical  root  finding  scheme.  However,  the  graphical  representation  of 
eq.  16a  by  Chang  [1]  helps  depict  the  propagating  and  evanescent  modes. 
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The  report  [l]'s  figure  number  4  is  repeated  here  as  Figure  3  for  convenier.  ;e . 
Equation  (16a)  is  rewritten  as  F^(k;a) 


k i  =  F„(k;a)  ;  F  (k;a)  =  ^-x(a) 
P  P 


(16b) 


Figure  3  shows  the  function  F  (ka)  for  p  *  0,  1,  2,...  as  a  function  of  a 

P 

for  a  dielectric  slab  of  premittivity  t  =  10eQ  and  thickness  d  of  1.27  mm  and 
desired  resonant  frequency  8GHz.  a  =  a(p)  is  then  given  by  intersecting  k£ 
with  the  set  of  curves  F^(k,a).  Only  a  finite  number  of  roots  appear  on  the 
real  a  axis.  In  this  case  for  p  =  0,  1,  2  we  have  a  =  3.05,  2.95,  1.075.  The 
rest  appear  as  a  continuous  mode  spectrum  of  the  evanescent  modes.  In  contrast, 
for  the  case  of  perfect  magnetic  walls  (perfect  cavity  -  no  radiation),  there 
are  not  only  a  fininte  number  on  the  real  axis  corresponding  to  propagating 
modes  but  also  eq.  (16b)  gives  an  infinite  set  of  discrete  roots  on  the 
imaginary  a  axis  corresponding  to  discrete  evanescent  modes  which  represent 
the  reactive  part  of  the  probe  impedance. 


11 


I 


Figure  3.  Graphical  Method  for  determining  the 

propagating  modes.  (After  Chang,  D.C.  [1]). 


2 


In  addition,  for  the  leaky  radiating  cavity  the  function  x(ct)  has  branch  points 

at  a  =  ±1  and  the  path  of  integration  is  deformed  upwards  (for  x  >  x  )  as 

o 

shown  in  Figure  4a.  The  integration  path  then  wraps  around  the  branch  cut  as 
shown  in  Figure  4b.  The  infinite  integral  in  eq.  (12),  would  reduce  to  the  sum 
of  the  branch  cut  integration  and  the  sum  of  the  residues  at  the  poles  on  the 
real  axis. 

Equation  (12)  can  be  written  as  follows 


13 


T 

Here  the  width  l  of  the  patch  is  chosen  such  that  only  the  dominant  TM^ 
mode  propagates. 

In  eq.  (20), 


X’(a)  = 


2kd 


(Mt(Y-D  *  2n\(-6c)  -  n2m2n)j_^ 


7tVn2-a2 

\»c  (Y-l)  ♦  2Q0(«m)  -  2»(2>.)}  (20  *  (^yj 

(  a  .  a  InQjq2-!  kd)  \ 

^r  | (cr2-l)  (n2-a2)  } 


-  Mr  |202„(ra,  *  ^ 


(21) 


The  physical  interpretation  of  the  poles  and  branch  cuts  is  as  follows.  The 
poles  on  the  real  axis  correspond  to  the  propagating  modes  which  define  the 
resonance  for  the  patch.  For  a  wide  patch  there  may  be  more  than  one  such 
guiding  mode.  In  the  case  of  perfect  magnetic  walls  a  series  of  poles  occur 
on  the  imaginary  axis.  These  correspond  to  the  discrete  evanescent  modes  ..a 
the  cavity  which  add  up  to  give  the  reactive  part  of  the  probe  impedance 
(see  Figure  5a).  However,  for  the  case  of  leaky  magnetic  walls  when  the  re¬ 
flection  coefficient  is  no  longer  unity  and  is  now  given  by  eq.  (14),  the 
evanescent  modes  appear  as  the  continuous  mode  spectrum  in  the  form  of  the 
branch  cut  integration  from  i®  to  o  (see  Figure  5b).  If  T  =  1  in  this  region, 
then  the  discrete  evanescent  modes  predicted  by  the  cavity  model  can  be  obtained. 

When  x  ^  0.  «  is  obtained  from 
P 

^  n2  +  -  x(iap)  =  P7*  (22) 

and  when  X  =  0»  that  is  T  =  1,  a  /  is  obtained  from 

P 

^  n2  +  a  2  =  pn  (23) 

P 
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It  can  be  shown  analytically  that  for  d/i  <<  1  the  difference  in  the 
above  two  cases  is  immaterial  for  integration  along  the  imaginary  axis 
and  we  can  then  put  T  =  1  in  this  region. 

Further,  the  modes  which  are  guided  down  the  patch  leak  out  as  spurious 
radiation  at  each  bounce  in  the  x  direction.  This  radiation  is  represented 
by  the  branch  cut  integration  from  0  to  1.  The  integrand  f(a)  given  in 
eq.  (18)  is  purely  reactive.  It  can  be  shown  analytically  that  this  spurious 
radiation  is  of  the  order  of  kd  which  can  be  ignored  since  the  resonant 
term  is  of  the  order  of  1/kd  at  resonance. 

Hence  the  field  under  the  infinite  microstrip  patch  antenna  reduces  to  the 
residue  at  the  pole  corresponding  to  the  propagating  mode  p=0  and  the  sums 
of  the  residues  at  the  poles  on  the  imaginary  axis  corresponding  to  the 
evanescent  modes  p  =  1,2,... 


3.0  TRUNCATED  MICROSTRIP  PATCH  ANTENNA 

Having  established  the  entire  spectrum  of  plane  waves  emanating  from  the 
exciting  probe  for  an  infinite  patch,  the  patch  can  now  be  truncated  to  bounce 
these  guiding  and  evanescent  modes  in  the  x  direction.  The  geometry  of  the 
truncated  patch  is  given  in  Figure  6.  The  length  of  the  patch  is  h  and  the 
probe  is  located  at  (x  ,  y  ) . 


Integration  path 


Figure  4a.  Deforming  the  path  of  integration. 


Propagating  and  evanescent  modes  for  perfect  magnetic 
walls . 


Leaky  cavity  with  an  angularly  dependent  reflection 
coefficient . 


Bouncing  the  mode  in  the  x  direction  would  amount  to  expanding  the  term 
ika(x-x  ) 

e  o  in  eq.  (18).  This  is  equivalent  to  the  sum  of  the  incident  waves 

and  the  waves  travelling  in  the  ±x  directions,  that  is 


:ikalx-x0l  +  Ceika(!  '  X)  +  Deik“(X  +  l) 


As  before  the  amplitudes  C  and  D  can  be  determined  at  h/2  and  -h/2.  In 
order  to  account  for  the  bounces  from  the  two  end  walls  of  a  microstrip  we 
need  to  replace  exp  (ika|x-XQ|)  in  eq.  (19)  and  (20)  by 


ika  x-x 
e  oi  + 


♦  [cos(ka(x+x  )) 

\-rH4^Z)el2kah  1  0 

+  f (Vn^-a^l  eikak  cos(ka(x-xQ))J 


^ikajx  xq|  eq#  £25)  is  replaced  by  eq.  (26)  to  obtain  the  entire  field 
under  the  truncated  microstrip  patch  antenna. 


4.0  INPUT  IMPEDANCE  TO  PROBE-FED  MICROSTRIP  PATCH  ANTENNAS 

The  electric  field  at  the  feed  probe  is  effectively  that  averaged  over  the 
area  of  the  probe,  i.e., 

1  2v 

<E  >  =  — —  /  dli  E  (x,y)  x  =  x  +  acosij),  y  =  y  +  asin<}>  (27) 

z  a  2it  '  z  o  o 

o 


where  a  is  the  probe  radius.  Utilizing  eqs.  (24)  and  (26), 


<E  >  = 

z  a 


W  -  a2  x'(a)  +  kid  I  l-r2(/n2_a2  ) 


2r(/n2  -  a2  >e 


i2kah 


|cos(2kaxQ)  +T(/n2  _  a2  )eka^ 

4tt1Aq  ®  cos  (2kaxn)  +  e^ 
ki  0  sin(kah) 

D  — Z 


P=2 

(even)  L 


®-*2  (Vr) 
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As  before,  we  consider  an  infinitely  long  microstrip  patch  and  establish 
all  the  propagating  and  evanescent  modes.  Substituting  T(a)  =  1  in  eq.  (12). 


Simplifying 


E 

z 


2jSo  eital1'’1. 
U  1  2a  e1?" 


i£7t  y 
£ 


0 


e 


i£7T 

£ 


+  e 


-ipn 
£ 


-ip7i  y /  -ipn  y 

£  £ 
e  ye 


+  e 


ipn 


i£n  y 
£ 


e 


a  =  a (p) 


(35) 


This  can  also  be  obtained  by  setting  T(a)  *  1  in  eq.  (24). 
For  p  even 


_  4nAo 
z  -  k£ 


2 

p  even 


cos  y )  cos  (&  y^j 


(36) 


For  p  odd 


-  2  e^k<*  x  xj  Sj_n  / pny  \  sin  ( pn  y  \ 

z  "  kl  p  odd  5  V  £  /  \^  £  7 


(37) 


These  eqs.  (36)  and  (37)  give  the  propagating  and  evanescent  modes  as  the  plane 
waves  bounce  in  the  y  direction  in  an  infinitely  extending  patch  with  perfect 

•  ika i v— x  i 

magnetic  walls.  Truncating  the  patch  and  expanding  e  |  where 

r(>^n^  -  a2)  -  1 

^ikajx  xQ|  then  replaced  by 


;iko,|x‘xoj  + 


ikah 


1-e 


i2kah 


jcos  ^ka(x+xQ)^ 


_  ikah 
+  e  cos 


(ka(x-xQ) 


) 


(38) 


kVnz-a7£=p n 
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The  sum  in  eq.  (38)  will  be  replaced  by  the  function  $(<*)• 


This  function  <j>(a)  should  satisfy  the  following  equations;  the  wave 

^^2  +  k2  az^j  0  =  i2ak6  (x-xq) 

and  the  boundary  condition 

4>’  (or)  =  0  at  x  =  ±h/2 

Assuming  a  solution  of  the  form 


♦  (<*)  =  2  An 


for  n  =  1, 


2... 


which  satisfies  the  orthogonality  relation 

h/2  Amcos  (k/2  X)  1  VOS(e7z  X)  dx  =  constaat  for  m  =  n 

-h/2  V 


substituting  (41)  into  (42)  we  have 


and 


equation 

(39) 

(40) 

(41) 

(42) 


I 


This  gives. 


♦(«>  =  i 


(im  \  (2m  \ 

co%~  \)  cos{—  xj 


(44) 


Here  a  satisfies  the  relation  k yjn^-a'^Z  -  p/i 
and  for  p  even,  where  p  =  2p' 

k2  «2  =  k2Q2 


Therefore, 

<H(a) 


(45) 


i4ak  £  _ 

h  m  k2n2 


(2mn  x  \  i 

( 2mnx \ 

r?  cos | 

— r —  oj  cosi 
V  h  /  ’ 

k  h  ) 

(46) 


and,  for  p  even 


E  =  ~ 


il6nAo  T  _ 

cos  | 

Vs  y’ 

LlJ 

1 

\  ( 2mn  x 

1  coshr 

^  cosi 

(tx)’ 

£h  ,  Z 

P  n 

k‘  .. .  (2£^* . 

Z - 

(47) 


Averaging  the  input  field  over  the  area  of  the  probe  we  have, 

fe  y°)  -  (¥  -o) 


/p  \  —  1 16/iAo  ^ 

<zz\  -  T  z,  l 

p  m 


cos  | 


2  n 
o 


/  cos( (x0  +  3  cos  $)^cos^-p  (yo  +  a  sin  $))  d<(i 


(48) 
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so  that  the  input  impedance  is  given  by 


The  modes  corresponding  to  p'  and  m  are  the  same  as  that  obtained  in  the 
cavity  model.  The  Bessel  function  is  the  converging  factor  similar  to  G 
in  [6]. 


6.0  RESULTS 


The  input  impedance  for  a  rectangular  patch  of  dimension  l  =  6.858  cm 
and  h  =  4.14  cm  mounted  on  a  grounded  dielectric  of  thickness  d  =  0.1588  cm 
and  er  =  2.5  and  excited  by  a  probe  of  diameter  2a  =  0.132  cms  has  been 
obtained  using  eqs.  (29)  and  (30a).  This  is  shown  in  Figure  7  as  a  function 
of  frequency.  The  resonant  resistance  has  been  compared  with  that  obtained 
by  Carver  and  Coffey  [10].  There  is  a  1%  difference  in  the  resonant  frequencies 
obtained  by  the  two  theories.  This  can  be  explained  as  follows. 


From  the  form  of  the  reflection  coefficient  T(a)  one  can  define  the  apparent 
end  admittance  of  the  microstrip  structure  as 


Y (k, a)  =  Y 

o 


i  -  r(q) 
l  -  r(a) 


where  Y  =  l -  1  is  the  characteristic  admittance  of  the  TEM  wave 

o  \  n  / 

in  a  parallel-plate  waveguide.  It  suffices  to  note  that  the  end  admittance 
in  general  is  a  function  of  both  frequency  and  angle  of  incidence.  Never¬ 
theless,  for  thin  slabs,  it  can  be  approximated  by  the  value  at  normal  inci¬ 
dence  at  the  end  walls  x  =  0  and  h.  Hence  the  concept  of  a  constant  end 
admittance  independent  of  the  angle  of  incidence  as  in  the  cavity  model  has 
merit  in  this  situation  provided  one  recognizes  the  dynamic  nature  of  the 
wall  admittance. 
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Now,  since  the  magnitude  of  the  reflected  wave  is  very  close  to  unity  we 
can  show  that 


G  =  (gYo)*  =  v(i)  where  v  = 

The  real  part  g  of  the  normalized  admittance  varies  linearly  with  kd  so  that 

g  is  now  independent  of  the  frequency.  When  g  is  plotted  as  a  function  of 

kd  shown  in  Figure  8  (repeated  from  [1]),  v  is  determined  by  the  slope  of  g 
-3 

and  is  8  x  10  .  On  the  other  hand,  the  conventional  theory  which  assumes 

that  the  end  conductance  can  be  calculated  approximately  from  an  open  slot 

_3 

radiation  in  the  absence  of  the  substrate  gives  a  slope  of  8.33  X  10  .  The 

error  in  omitting  the  surface  wave  radiation  is  typically  of  the  order  given 
by  (1  -  e^^kd  and  is  not  significant  when  the  substrate  is  thin. 

As  for  the  end  susceptance  of  a  microstrip  of  width  we  can  first  define 
an  equivalent  length  AS  of  an  open  circuit  according  to 

B  =  (bY  )S  =  uiz  e  AS  where  =(“rr)  t 
o  r  o  d  Vnkd/  d 

which  is  also  plotted  in  Figure  8.  Unlike  the  quasi-static  theory  which 
gives  a  result  independent  of  frequency,  the  dynamic  nature  of  our  result 
is  very  explicit.  Experimental  results  by  Johnk  and  Chang  [1]  show  that 
the  Wiener-Hopf  theory  properly  incorporates  the  high  frequency  effects 
into  the  modeling  of  the  open  end  susceptance  whereas  the  static  approach 
is  inadequate.  An  especially  encouraging  result  of  this  is  the  fact  that 
the  incorporation  of  the  angle  of  incidence  in  the  Wiener-Hopf  reflection 
coefficient  produces  a  susceptance  that  is  in  excellent  agreement  with  the 
experimental  results. 

In  view  of  the  above  discussion  it  is  obvious  that  the  discrepancy  between 
the  two  results  shown  in  Figure  7  arises  due  to  the  difference  the  wall 
susceptance  used  by  the  two  theoretical  models.  Further,  the  dynamic  nature 
of  the  wall  susceptance  makes  the  present  theory  more  appropriate  for 
millimeter  wavelength  antennas. 
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Figures  9  and  10  show  the  resonant  resistance  and  reactance  for  the  dominant 
mode  as  the  substrate  thickness  varies.  As  expected  the  resonant  frequency 
shifts  and  the  beamwidth  increases  with  increasing  electrical  thickness. 
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kd  (NORMALIZED  SUBSTRATE  THICKNESS) 

Plot  of  Real  Part  of  Normalized  Admit!. ance 
vs.  Normalized  Substrate.  Thickness. 

(After  Chang,  D.C.  [1]). 


7.0  CONCLUSION 


This  report  has  presented  a  general  theory  for  the  analysis  of  the  input 
impedance  to  a  probe-fed  microstrip  patch  antenna  where  the  entire  spectrum 
of  plane  waves  emanating  from  the  exciting  probe  has  been  considered  and 
an  angularly  dependent  wall  reflection  coefficient  has  been  incorporated. 

The  primary  emphasis  to  its  application  is  for  the  case  of  millimeter  antennas 
where  the  substrate  is  electrically  thick  and  makes  it  necessary  to  include 
the  dispersion  effects  of  the  dielectric. 

It  is  seen  here  that  there  are  two  kinds  of  waves  which  originate  from  the 
probe;  those  which  are  evanescent  and  decay  in  the  vicinity  of  the  probe 
therefore  constituting  the  reactive  part  of  the  impedance  and  those  which 
propagate  down  the  patch  and  are  the  guiding  modes.  The  guiding  modes  are 
discrete  and  finite  in  number  tor  a  given  patch  size  and  desired  frequency 
of  operations  whereas  the  evanescent  modes  form  a  continuous  spectrum. 

Further,  the  magnitude  of  the  spurios  radiation  which  leaks  out  of  the 
cavity  as  the  entire  spectrum  of  plane  waves  bounce  back  and  forth  has  been 
analytically  estimated  to  be  of  the  order  of  kd.  The  continuous  spectrum  of 
evanescent  modes  for  the  thick  substrate  or  leaky  walls  is  in  contrast  to  the 
case  of  the  perfect  magnetic  walls  for  which  it  can  be  shown  analytically  and 
confirmed  by  preliminary  numerical  analysis  that  the  evanescent  modes,  like 
the  propagating  modes,  are  discrete.  Hence  the  total  field  reduces  to  merely 
calculating  the  residues  at  the  poles  corresponding  to  all  the  modes  since 
there  is  no  (by  assumption)  radiation  from  the  cavity.  Although  this  assump¬ 
tion  of  discrete  modes  is  also  reasonable  for  substrates  which  are  electrically 
thin,  the  radiation  has  to  be  accounted  for  in  substrates  which  are  electrically 
thick  as  in  the  case  of  millimeter  wave  microstrip  antennas.  Analytically  it 
is  shown  that  for  the  limiting  case  of  perfect  magnetic  walls  this  general 
theory  reduces  to  a  form  similar  to  that  obtained  in  the  cavity  model.  In  order 
to  do  this,  the  single  summation  in  eqs.  (34)  and  (35),  which  includes  the 
entire  spectrum,  was  made  to  take  the  form  of  the  familiar  double  summation  which 
appears  in  the  cavity  model  analysis.  It  should  therefore  be  emphasized  that 
establishing  the  mode  spectrum  for  an  infinitely  long  patch  and  then  bouncing 
these  modes  for  the  truncated  patch  not  only  gives  all  the  modes  in  both 
directions  under  the  patch  but  also  makes  the  numerical  and  mathematical  analysis 
simpler. 
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